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We demonstrate experimentally the usefulness of selec- 
tive pulses in NMR to perform quantum computation. Three 
different techniques based on selective pulse excitations have 
been proposed to prepare a spin system in a pseudo-pure 
state. We describe the design of novel "portmanteau" gates 
using the selective manipulation of level populations. A selec- 
tive pulse implementation of the Deutsch-Jozsa algorithm for 
a two-qubit and a three-qubit quantum computer is demon- 
strated. 

PACS:03.67.Lx,76.70.-k 



I. INTRODUCTION 

The idea of exploiting the intrinsically quantum me- 
chanical nature of physical systems to perform compu- 
tations has generated a lot of excitement recently [|l,||. 
Logical operations in quantum computation are imple- 
mented on quantum bits(qubits), where a qubit can be 
any quantum two-level system. The two eigenstates are 
mapped onto logical and 1. While all classical compu- 
tation can be performed using the above mapping, the 
fact that a qubit can exist in a general coherent superpo- 
sition of logical states and 1, leads to new possibilities 
for computation. The realisation that computation can 
be performed reversibly oH , paved the way for the quan- 
tum mechanical implementation of logic gates through 
unitary transformations P,p|. The power of quantum 
computing lies in the fact that a single input state of a 
quantum computer can be a coherent superposition of all 
possible classical inputs. Consequently, algorithms that 
are intrinsically quantum in nature can be designed, to 
solve problems hitherto deemed intractable on classical 
computers R-O]. A major hurdle in achieving quantum 
computing experimentally, is that of preserving quantum 
coherence while the computation is being performed, and 
the search for such ideal quantum systems yielded single 
charged ions confined in an ion trap pO| | and nuclear spins 
in a liquid as possible choices. 

It has been demonstrated that assemblies of nuclear 
spins in a liquid, which are largely isolated from their en- 
vironment and have long relaxation times (so coherence 



is retained for a while), can be used to build quantum 
information processors m . A system of N spins can ex- 
ist in entangled quantum superposition states and can 
be thought of as an TV bit quantum computer. How- 
ever, quantum computing requires pure states as inputs, 
whereas nuclear spins at thermal equilibrium are in a 
statistical mixture of pure states. It was demonstrated 
recently that it is possible to perform quantum comput- 
ing with mixed state ensembles rather than on an iso- 
lated system in a pure state |11[ . The problem is cir- 
cumvented by creating within the overall density matrix 
of the system, a sub-ensemble that behaves like a pure 
state. Techniques to prepare such "pseudo-pure" states 
have been proposed by different groups O-IIsI . Previous 
workers in the field have employed various NMR methods 
like non-selective pulses, rf gradients, coherence transfer 
via J-coupling and simultaneous multi-site excitation to 
create pseudo pure states, construct universal logic gates 
and implement quantum algorithms for two and three- 
qubit systems |llj-|ll. 

In this paper, we explore the utility of transition- 
selective and spin-selective pulses, and exploit the non- 
commuting nature of operations on connected transitions 
to prepare a spin system in a pseudo-pure state, exe- 
cute different logical operations simultaneously and im- 
plement the Deutsch-Jozsa quantum algorithm on a two 
and a three-qubit system. The Ti relaxation times in the 
molecules used, is of the order of a few seconds (3.4 sees - 
4.6 sees), whereas T2 relaxation occurs within an interval 
of around 1 sec. Selective excitation has been achieved 
using low power, long duration pulses of a rectangular 
shape. The length of these pulses is tailored to achieve 
sufficient selectivity in the frequency domain without per- 
turbing the nearest line, and hence depends on the mag- 
nitude of the smallest J coupling present. The duration 
of the pulses applied varies from 100 ms to 263 ms (for J 
couphngs of 9.55 Hz to 3.8 Hz). For small computations, 
such as the ones performed here, drastic decoherence or 
dephasing does not occur during the duration of these 
selective pulses. However, the deleterious effects of such 
selective pulses must be considered and compensated for, 
whenever larger computations are attempted. 



II. CREATION OF PSEUDO-PURE STATES 

The logical labeling technique to create pseudo-pure 
states is broadly categorised by the fact that unitary 
transformations are used to redistribute the populations 
of states, such that an effective pure state is obtained in 
the sub-manifold of qubits(spins) to be used for computa- 
tion, and ancillary qubits are used as "labels" . While the 
concept underlying the logical labeling method of pseudo- 
pure state creation has been delineated by Chuang et. 
al. ||l3|,|l6|, there have been very few experimental im- 
plementations of such an elegant technique |2^. We 
have designed a few novel pulse schemes using transition- 
selective pulses to create such logically labeled pseudo- 
pure states. Consider a three spin-1/2 system(AMX), 
with the energy levels labeled as in Figure IT]. 



homonuclear three-spin system of 2,3 dibromopropionic 
acid is shown in Figure |^. 
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FIG. 1. The creation of a pseudo-pure state in an AMX 
three spin system using logical labeling, (a) The population 
distribution of the thermal equilibrium state, (b) The popu- 
lation distribution of a pseudo-pure state, created by invert- 
ing the populations of the two single-quantum A transitions 
shown in (a) by long arrows. 

The selective inversion of two unconnected single- 
quantum transitions of the A spin(| tti) ^ I iti) and 
I Tit) — * I iit)) would lead to the creation of a logi- 
cally labeled pseudo-pure state, with A being the "label 
qubit" and M, X being the "work qubits" available for 
computation. The first four eigenstates (labeled by the 
first spin being in the | t) state) now form a manifold 
that corresponds to a two-qubit pseudo-pure state while 
the last four(labeled by the first spin being in the | |) 
state) form a separate manifold that corresponds to an- 
other two-qubit pseudo-pure state. The creation of a 
pseudo-pure state by this method leads to relative pop- 
ulation differences of 
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FIG. 2. A logically labeled pseudo-pure state in the 
homonuclear three spin system(AMX) of 2,3 dibromopropi- 
onic acid, (a) The equilibrium proton spectrum is shown, with 
the three protons labeled A, M and X resonating at 5a = 3.91 
ppm, Sm = 3.69 ppm and Sx = 4.48 ppm respectively, (b) 
The selective inversion of the two (nearly overlapping) central 
transitions of the A spin leads to the creation of a logically 
labeled pseudo-pure state, which has been read by a small 
flip angle(10'') detection pulse. Long, low-power rectangular 
pulses have been used for selective excitation. 

The pseudo-pure state has been distilled by manipulat- 
ing unconnected single quantum transitions of the label 
qubit A, as detailed in Equation (2.1). The transition- 



Homonuclear three spin case: The experimental 
creation of a logically labeled pseudo-pure state in the 



selective tt pulses were applied on the two central (nearly 
overlapping) transitions of the A spin. 
Heteronuclear three spin case: The experimental 
creation of a logically labeled pseudo-pure state in the 
heteronuclear three spin system of 4-fluoro,7-nitro ben- 
zofurazan is shown in Figure g[ Two selective tt pulses 
were applied on the central, nearly overlapping uncon- 
nected single-quantum transitions (| TTT) ~^ I TTT) £^nd 
I TTT) — * I TiT)) of the A spin(the proton in this case). 
The A spin is the "label" qubit and the other two spins 
(the third spin being fluorine in this case) are the "work" 
qubits. It is to be noted that the spectral pattern of the X 
spin in the two pseudo-pure states created in the homonu- 
clear and the heteronuclear systems are mirror images of 
each other i.e. intensities of (0,2,2,0) are obtained for 
the X spin in Figure 0(b) while the X spin multiplet pat- 
tern is (2,0,0,2) in Figure |3[ This difference reflects the 
relative sign of the coupling constants in these systems. 
Heteronuclear ^^F-^H spin systems are useful for quan- 
tum computing as they have the twin advantages of good 
sensitivity and long relaxation times. 

Other methods to implement a logically labeled 
pseudo-pure state can be designed, based on the selective 
manipulation of the populations of multiple-quanta. For 
instance, the inversion of the double quantum (| Ttt) -^ 
I TTT)) I followed by the inversion of the single quantum 



transition (| fjl) — > | Hi)), leads to another pseudo- 
pure state. The redistribution of equiHbrium populations 
leads to relative population differences for the pseudo- 
pure state 
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Experimentally, the zero quantum can be inverted by a 
cascade of tt pulses on two regressively connected single 
quantum transitions by the cascade tt^'^ tt^-^ tt^'^ (spectra 
not shown) [ P2[ . 

The state of the system after the creation of the log- 
ically labeled pseudo-pure state has been read out by a 
small flip angle detection pulse in each case. While it is 
usual in quantum computing to use pulses of flip angle 
90° for the read out operation, this will not provide a 
useful output for logical labeling experiments. To illus- 
trate this point, consider a logically labeled pseudo-pure 
state for the three spin AMX system, with A being the 
label qubit, and M and X the work qubits (Eqn. [2.1[ ). 
The traceless, deviation density matrix corresponding to 
this pseudo-pure state can be described in terms of prod- 
uct operators (which are just products of spin angular 
momentum operators) 



^ p~pure 



= M, +X, + 4A,M,X, 



(2.4) 



A detection pulse of flip angle a leads to the NMR ob- 
servable terms 

(M^ + X^) sina + {4A^M,X, + 4A,M^X, + 

iA^M^X^) cos^ asina (2.5) 

Hence, a 90° detection pulse would not be able to read 
out all the product operators present in the density ma- 
trix (Eqn. ^.4|), and a small angle read pulse is required. 



III. NOVEL QUANTUM LOGIC GATES 



FIG. 3. The creation of a logically labeled pseudo-pure 
state in 4-fluoro,7-nitro benzofurazan(AMX system). Tran- 
sition selective n pulses have been applied on the two over- 
lapping central transitions of the "label" qubit(the proton A). 
The asterisk labels the solvent peak. The normal proton and 
fluorine spectra are shown in (a) and (c) with the coupling 
constants Jam = 8.1Hz, J ax = 8.01Hz and Jmx = 3.81Hz 
and the spins resonating at Sa = 7.28ppm, Sm = 8.57 ppm 
and Sx ~ —110.35 ppm respectively. The proton and fluorine 
spectra corresponding to the pseudo-pure state are shown in 
(b) and (d) respectively, with the "work" qubits M and X 
being in a logically labeled pseudo-pure state. The state of 
the system has been monitored by a small flip angle detection 
pulse. 

Experimentally, the double quantum can be inverted 
by a cascade of tt pulses on progressively connected sin- 
gle quantum transitions |22| . The inversion of the sin- 
gle quantum transition | jtt) -^ \ iti)i followed by the 
inversion of the zero quantum | |J,|) -^ \ jti); leads 
to yet another pseudo-pure state. The redistribution of 
equilibrium populations leads to the relative population 
differences 



TTT TTi TiT Tii TTT iTi iiT iii 

3/2 1/2 1/2 1/2 -1/2 -1/2 -1/2 -3/2 



(2.3) 



We now proceed towards the implementation of quan- 
tum logic gates using NMR lO jll] • The two-qubit quan- 
tum XOR(or controlled-NOT) gate has been demon- 
strated to be fundamental for quantum computation 
and has been implemented in NMR by a selective [tt 
pulse on a single transition ||ll|] . It has been proved that 
the reversible quantum XOR gate, supplemented by a 
set of general onc-qubit quantum gates, is sufficient to 
perform any arbitrary quantum computation |q| . Hence, 
with the achievement of the experimental implementa- 
tion of the XOR and one-qubit gates in NMR, it does not 
seem necessary to look for the design and construction of 
other gates. Nevertheless, it is important in terms of 
complexity in large circuits and ease of experimental im- 
plementation, to look towards the design of efficient logic 
networks. In this direction, gates that achieve the imple- 
mentation of two or more logic operations simultaneously 
would be useful in reducing computational time in cir- 
cuits that require a large number of logical operations. 
We detail the design and experimental implementation 
of such gates here; borrowing from Lewis Carroll, we call 
such "many-in-one" gates "Portmanteau" gates p^ . 
The logical SWAP operation: Consider a two-spin 
system (AX) with each spin being a qubit, the spin A 
being the first qubit and the spin X, the second qubit. 



The eigen states of this system can be represented by 
|ei,e2), where ei and £2 are or 1. The logical SWAP 
operation completely exchanges the states of a pair of 
qubits, from |ei, £2) to |e2, ei), the unitary transformation 
corresponding to such an operation being 



[/«■ 
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10 
10 
1 



(3.1) 



This gate might be useful during the course of a com- 
putation when qubits need to be permuted pq |. In spin 
systems where some scalar J couplings are ill resolved, 
the logical SWAP could be used to compensate for the 
missing couplings pO(| . 

Madi et. al. pffl, have implemented the SWAP oper- 
ation using an INEPT-type sequence, with non-selective 
rf pulses and J-evolution. It is interesting to note that 
the logical SWAP operation can be achieved by selec- 
tively interchanging the populations of the zero quantum 
levels. Since these levels are not connected by single- 
quantum transitions, the population exchange will have 
to be achieved indirectly. 

The inversion of the zero quantum (SWAP) requires 
a cascade of three selective tt pulses on regressively con- 
nected transitions, for example, [7r];^^[7r]^i [Tr];^^ P3]. 



tions of the XOR, XNOR and NOT gates (all with their 
outputs on the first qubit) are given by 



XOR 



ei, £2) — > |ei ® £2, £2); t^xoR = 



I > XNOR I — - — > „i 

|ei,e2) — > |ei©e2,e2); U^^ 



I > NOT |_ , 1 



10 
1 
10 
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10 
10 
10 
1 

10 
1 
10 
10 



(3.2) 



The superscript 1 indicates that the output of the gate 
is obtained on the first qubit. The matrices C/xor ^t*^- 
corresponding to the output on the second qubit can be 
similarly constructed. 

Logical SWAP+XOR(XNOR): The execution of a 
logical SWAP operation followed by an XOR gate (with 
its output on the first qubit), can be defined through 
its action on |ei,e2), and leads to the final state |ei © 
£2,61) (or to |ei © e2,ei} for a SWAP followed by an 
XNOR gate with its output on the first qubit) with their 
explicit matrix representations being 




7.6 7.4 7.2 7.0 6.8 6.6 [ppm] 
FIG. 4. The logical SWAP operation implemented on the 
two-spin system of Coumarin. The result of the application 
of a [■k]x pulse on one of the X transitions, followed by a [n]x 
pulse on the regressively connected A transition is shown, 
which corresponds to a logical SWAP on the non-equilibrium 
state of the spin system. A small angle (lO'^) read pulse has 
been used. The ID spectrum and structure of Coumarin is 
shown in Figure n. 

Since this would lead to the same spectrum as the 
equilibrium spectrum, a non-equilibrium state has been 
first created by preceding the cascade with a selective 
[tt]^! pulse, yielding the cascade H^^^lf ' Wi' Wx' = 
[7r];^4^]^^- This amounts to the execution of a logical 
SWAP operation on a non-equilibrium state. The ex- 
perimental implementation of this operation is shown in 
Figure on the two spin system of Coumarin. 

We now explore the implementation of gates that re- 
alise various combinations of the SWAP, XOR(XNOR) 
and NOT operations. The action and matrix representa- 



U. 



SWAP + XOR. 



u. 



SWAP-fXNOR 



^XOR *-^SWAP 



t-^YWOR ^S' 
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1 



(3.3) 



An implementation of these operations requires the ap- 
plication of selective ['k\x pulses consecutively on two re- 
gressively connected transitions, and the resulting spec- 
trum on two bits is identical to Figure ^ 
XOR(XNOR)+Logical SWAP+NOT: The imple- 
mentation of an XOR gate (with the output on the first 
qubit), followed by a SWAP operation and then a NOT 
gate on the first qubit, corresponding to a final state of 
|e2, ci ©£2) can be experimentally achieved by transition- 
selective [7r]a; pulses applied consecutively on two progres- 
sively connected transitions. Reversing the order of ap- 
plication of the pulses leads to the final state |ei © e2, eT) 
which corresponds to an XNOR-|-Logical SWAP-I-NOT 
gate, with the output on the second qubit. 

These gates correspond to the unitary matrices given 
by 



U: 



XOR+SWAP+NOT 



■-^NOT '-'SWAP '-'XOR 



JJ TJ^ 

■ '-'.SWAP '-^ vr 



1 
10 
10 
10 



U: 



XNOR+SWAP+NOT 



JJ^ JJ 






10 
10 
1 
10 



(3.4) 

The experimental implementation of such gates is 
shown in Figure ||(b),(c). It is interesting to note that 
these operations do not commute, so the order in which 
the pulses are applied is important and its reversal leads 
to different logical operations. 

NOT+Logical SWAP: The NOT gate followed by a 
logical SWAP operation on two qubits (or vice versa since 
these operations commute) leads to | £2 , ef) when applied 
to the state |ei, £2). This action suffices to determine the 
unitary matrix for the above operation, given by 



c/, 



NOT+SWAP 



f^SWAP ^NOT "~ f^NOT ^ S^ 



1 
10 
10 
10 



(3.5) 

The experimental implementation has been achieved 
by selectively inverting the populations of the double- 
quantum levels. A cascade of transition-selective [tt]^ 
pulses has been applied on two progressively connected 
transitions (Figure ||(d)). 

The implementation of various portmanteau gates on a 
thermal initial state is shown in Figure ra for the two-spin 
system of Coumarin. The same pulse schemes are ex- 
pected to implement the above logic operations on other 
initial states (for instance, a pseudo-pure or a coher- 
ent superposition of states) as well. As an illustration, 
consider the portmanteau gate XOR-hSWAP+NOT. The 
unitary matrix corresponding to it (Eqn. 3.4) can be de- 
composed into two matrices 



a 



ui 



XOR+SWAP-l-NOT 

10 

10 

i 

i 



^u^^u^^- 
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i 
1 



(3.6) 



The matrices U:^} and U^^ correspond to selective ['k]^ 
pulses on the Ai and X2 transitions respectively. A 
low-power, long duration rectangular pulse is applied to 
achieve the desired selectivity. These transition-selective 
pulses can be expanded in terms of single-transition op- 
erators (the expansion is independent of the state of the 
spin system P2| , p4| ). This realisation in terms of sin- 
gle transition operators is valid when the power of the rf 
pulse is low compared to the the J-coupling and the chem- 
ical shift difference between the spins {wi « 2tiJ « 
5 ax)- One is thus able to realise the unitary transfor- 
mations required to implement the desired logical opera- 
tions, without prior knowledge of the state of the system. 



However as noted recently [g5[, a more complete Hamil- 
tonian description might be required to describe evolu- 
tion under selective pulses, in order to fully establish the 
generality of the above schemes. 
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FIG. 5. The experimental implementation of various port- 
manteau gates, (a) Reference spectrum of Coumarin at ther- 
mal equilibrium, (b) The implementation of an XOR -|- 
SWAP + NOT (by the pulse cascade W^i followed by [7r]f=, 
where Ai and X2 refer to progressively connected transitions 
of spins A and X respectively) . (c) The implementation of an 
XNOR + SWAP + NOT ([ttJ^^ followed by M^i). (d) The 
implementation of a NOT + SWAP ([7r]^i[7r]^2[7r]^i). The 
state of the spin system is read by a small (lO") angle pulse 
in each case. 

It is to be noted that we have used selective pulses 
polarised along the x direction. The phases of the pulses 
are important, and have been experimentally ensured by 
appropriate phase cycling schemes. 

IV. IMPLEMENTATION OF THE 

DEUTSCH-JOZSA(D-J) QUANTUM 

ALGORITHM 



Finally, we experimentally implement the D-J algo- 
rithm using selective pulses. The D-J algorithm deter- 
mines whether an unknown function f{x) is constant or 



balanced ||]. In the simplest version, f{x) is a mapping 
from a single bit to a single bit and the function is con- 
stant if f{x) is independent of x and it is balanced if f{x) 
is zero for one value of x and unity for the other value. 
The generalisation to N bits is conceptually simple and 
f{x) in this case is constant if it is independent of x and 
balanced if it is zero for half the values of x and unity 
for the other half. A classical computer proceeding deter- 
ministically would require up to 2^^^ + l function calls to 
check if f{x) is constant or balanced; even if half the in- 
puts have been evaluated and all outputs have been found 
0(or 1) one cannot conclude that the function is constant. 
The quantum version of the algorithm determines if the 
function is balanced or constant using only a single func- 
tion call. For the one bit case this is achieved by evalu- 
ating the value of /(O) © /(I) (where © denotes addition 
modulo 2). The binary function / is encoded in a unitary 
transformation by the propagator Uf by including an ex- 

tra input qubit such that \x)\y) — > \x)\y (B f{x)). The 
four possible functions for the single-bit D-J algorithm 
are categorised as 
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The unitary transformations corresponding to the four 
possible propagators Uf can be easily constructed: 
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(4.1) 



The algorithm requires one input spin and one work spin. 
Using the propagator Uf and appropriate input states, 
one can proceed with the implementation of the algo- 
rithm. Previous workers in the field [^ 17|, used a com- 
bination of spin-selective (7r/2) pulses and evolution un- 
der the scalar coupling J, to encode the D-J algorithm 
on a two-qubit quantum computer. 

We have implemented the D-J algorithm using spin- 
selective and transition-selective n pulses. The exper- 
iment begins with both qubits in a superposition of 
states, achieved by a non-selective [Tr/2]y pulse on both 
spins. After application of the propagators Ui, the first 
qubit (the "control" qubit) remains in the superposition 
state while the desired result (/(O) © /(I)) is encoded 
as the appearance or disappearance of the lines of the 
target qubit. The Ui transformation corresponds to the 
unity operation or "do nothing" , while the U2 transform 
is achieved by a spin-selective [k]^ pulse on the con- 
trol qubit. The C/3 and U4 transformations are imple- 



mented by selective n pulses on the | |t) — > | ||) and the 
I TT) ~^ I Ti) transitions respectively. 
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FIG. 6. A selective pulse implementation of the D-J quan- 
tum algorithm on a two-qubit system 5-nitro furaldehyde, at 
room temperature on a 400 MHz spectrometer. The results 
after applying the unitary transformations Ui, U2, U3 and U4 
on a coherent superposition are shown in (1), (2), (3) and (4) 
respectively. 

The implementation of this algorithm does not re- 
quire pure initial states as similar results can be read 
out from the spectrum if one starts with thermal initial 
states instead. In a single measurement, one can distin- 
guish between constant and balanced functions on the 
basis of the disappearance of the lines of the target qubit 
in the spectrum. These predictions are borne out by 
the experimental spectra in Figure 0. The phase of the 
transition-selective pulse (to implement U3 and C/4) has 
been stepped through {x,—x,y,—y) to suppress phase 
distortions and leads to the total suppression of the tar- 
get qubit lines and the retention of only one line of the 
control qubit. 

The algorithm to distinguish between the two cate- 
gories (constant or balanced) of two-bit binary functions 
is implemented on a three-qubit NMR computer, by eval- 

uating |a;)|y)|z) — > \x)\y)\z®f{x,y)). The eight possible 
(2 constant and 6 balanced) two-bit binary functions are 
categorised as 
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Previous researchers used shaped pulses generated by 
an rf waveform generator to implement the two-bit D-J 
algorithm using three qubits; the pulse waveforms were 



tailored to selectively excite two or more frequencies si- 
multaneously [[l9| . 

We describe here a selective pulse implementation of 
the D-J algorithm using simple rectangular pulses (Fig- 
ure ItI). 
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the transition selective pulses have been stepped through 
(x, —X, y, —y) as in the two-qubit case and a similar logic 
prevails in explaining the spectral pattern obtained. Un- 
like the previous implementation of the three-qubit D-J 
algorithm |19|, the phase cycling here achieves a com- 
plete suppression of the multiplets of both the target 
spins when the function is balanced (/s — fs). 

It has been demonstrated that selective pulse tech- 
niques in NMR are a powerful tool to build quantum 
information processors. The distilling of a pseudo-pure 
state from a thermal state, the simultaneous implementa- 
tion of different logical operations to save computational 
time and a two and three-qubit implementation of the D- 
J quantum algorithm has been experimentally achieved 
using such techniques. 

Acknowledgements 

The use of the AMX-400 high resolution FTNMR spec- 
trometer of the Sophisticated Instruments Facility, In- 
dian Institute of Science, Bangalore, funded by the De- 
partment of Science and Technology, New Delhi, is grate- 
fully acknowledged. 



14) 



(3) 





4.5 4.4 4.3 4.2 4.1 4.0 3.9 3.8 [ppm] 

FIG. 7. Selective pulse implementation of the D-J algo- 
rithm on the three-qubit system of 2,3 dibromopropionic acid. 
The two constant functions /i and /2 are shown in (1) and 
(2) respectively, while the unitary transforms corresponding 
to the balanced functions /a — fs are implemented in (3)- (8) 
respectively. 

These low power, long duration transition-selective 
pulses are applied consecutively, and do not require any 
special hardware for their application. The unitary trans- 
forms have been implemented on a coherent superposi- 
tion of all the three qubits, achieved by a non-selective 
7r/2 pulse on a thermal initial state. The two constant 
functions /i and /2 correspond to the unity operation and 
a spin-selective [ttJj, pulse on the multiplet of the control 
qubit, respectively. The unitary transformations encod- 
ing the six balanced functions /a — fs are implemented 
by selective pulses on the transitions of the control qubit, 
taken two at a time i.e. the pulses can be described 
by [7r,7r,0,0], [^,0,^,0], [7r,0,0,7r], [0,7r,7r,0], [0,^,0,^] 
and [0, 0, tt, tt] on the four transitions where denotes 
no pulse on that particular transition. The phases of 
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